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ABSTRACT 
> 

(^ ' We consider the acceleration of charged particles near ultra-relativistic shocks, 



with Lorentz factor Fg ::^ 1. We present simulations of the acceleration process 



l^ , and compare these with results from semi-analytical calculations. We show 

that the spectrum that results from acceleration near ultra-relativistic shocks 
, is a power law, N{E) (x -E^^, with a nearly universal value s « 2.2 — 2.3 for 

Oh' the slope of this power law. 

O I We confirm that the ultra-relativistic equivalent of Fermi acceleration at 

C/3 . a shock differs from its non-relativistic counterpart by the occurence of large 

a 

anisotropics in the distribution of the accelerated particles near the shock. 



In the rest frame of the upstream fluid, particles can only outrun the shock 



JH I when their direction of motion lies within a small loss cone of opening angle 

9c ~ F^^ around the shock normal. 

We also show that all physically plausible deflection or scattering mech- 
anisms can change the upstream flight direction of relativistic particles orig- 
inating from downstream by only a small amount: A9 ^ TfT^- This limits 
the energy change per shock crossing cycle to AE ^ E, except for the first 
cycle where particles originate upstream. In that case the upstream energy is 
boosted by a factor ~ Fg for those particles that are scattered back across the 
shock into the upstream region. 

Key words: cosmic rays - acceleration of particles - shock waves - gamma- 
ray bursts. 
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2 Achterberg et al. 

1 INTRODUCTION 

Diffusive slioclc acceleration (DSA), independently proposed by various authors in the con- 
text of the acceleration of Galactic cosmic rays (Krymskii, 1977; Axford, Leer & Skadron, 
1978; Bell, 1978 and Blandford & Ostriker, 1978), has become a paradigm for the pro- 
duction of energetic particles near strong shocks in a magnetized plasma. It is believed to 
be the mechanism responsible for the production of Galactic cosmic rays with an energy 
E < 10^^ eV near the Sedov- Taylor blast waves associated with supernova remnants, below 
the so-called 'knee' in the cosmic ray spectrum. 

Diffusive shock acceleration relies on repeated scattering of charged particles by magnetic 
irregularities (Alfven waves) to confine the particles for some time near the shocks. This lim- 
its the mean free path of these particles to values much less than the one derived for cosmic 
rays in the general interstellar medium, Aism ~ 0.2 pc for protons at -E ~ 5 GeV, but still 
significantly larger than the thickness of the shock. This may explain the good correspon- 
dence between the radio maps of some supernova remnants, which show the distribution of 
the GeV electrons and the magnetic field responsible for the synchrotron emissivity, and the 
X-ray observations which show the distribution of the hot T ~ 10^ K shocked interstellar 
gas (Achterberg, Blandford & Reynolds, 1994; Aschenbach & Leahy, 1999). Direct observa- 
tional evidence for particle acceleration seems to be provided by the observation of X-ray 
synchrotron emission from electrons with E ~ 100 TeV in SN1006 (Koyama et al., 1995). 

One of the main strengths of diffusive shock acceleration is that it predicts a power-law 
spectrum which is not far from that required for a theory of Galactic cosmic rays. The 
observed spectrum of these particles, Nohs{E) dE oc i5~^'^ dE, is thought to arise from the 
effects of energy-dependent propagation on a source spectrum in the range Ns{E) dE oc 
E''"^ dE with power-law index s in the range 2.1<s<2.4. A similar spectrum is inferred 
for the GeV electrons responsible for the synchrotron emission in such non-thermal radio 
sources as supernova remnants and the lobes of radio galaxies, where the slope of the power 
law lies in the range s ~ 2 — 2.7. 

The formation of a power-law spectrum is the result of the competition between the 
energy gain per shock crossing cycle, say from upstream to downstream and back, and the 
chance of escape from the shock per crossing cycle, Pesc- Shock acceleration therefore is a 
realization of the statistical acceleration process proposed by Fermi (1949). The slope of the 
resulting power-law distribution is given by 
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Particle acceleration 



ln(l/Prct; 



Here Prct = 1 — Vcsc is the mean probability per cycle that a particle returns to the shock 
and re-crosses into the upstream medium, E^/Ei is the ratio of final and initial energy in a 
cycle, and the angular brackets indicate the average value. 

For non-relativistic shocks one can use simple kinematics, assuming elastic scattering 
of the accelerated particles in the upstream and downstream fiuid rest frames, which leads 
to near-isotropic particle distributions. The calculation of the mean energy gain and return 
probability is straightforward, and the resulting slope s depends only on the compression 
ratio r = Pd/Pu at the shock (e.g. Bell, 1978): 

.^^^ (2) 

r — 1 

Here Pu (Pd) is the fiuid density just upstream (downstream) of the shock. For a strong 
non-relativistic shock in an ideal mon- atomic gas one has r = 4 and s = 2. The fact that 
the observed spectra are often steeper is usually explained as the effect of energy- dependent 
escape of accelerated particles from the source (or the galaxy as a whole), or the modification 
of the shock by the back-reaction of the accelerated particles: the gradient in the pressure 
of the accelerated particles slows down and heats the incoming fiuid, decreasing the shock 
compression ratio. 

Reviews of the theory of diffusive shock acceleration at non-relativistic shocks can be 
found in Drury (1983), Blandford & Eichler (1987), Jones & Ellison (1991), Achterberg 
(1993) and Kirk (1994). 

In this paper, we consider the process of Fermi-type shock acceleration in the limit 
Fs ^ 1. We briefly review the importance of such shocks in astrophysics in Section 2. The 
basic kinematic constraints are discussed in Section 3. There we will show that the energy 
gain per shock crossing is of order unity, (Ef/Ei) ^ 2, except for the flrst shock crossing 
(Gallant & Achterberg, 1999). Possible scattering mechanisms are discussed in Section 4, and 
the maximum energy that can be achieved with and without losses is discussed in Section 
5. We present numerical simulations of the acceleration process in Section 6, and compare 
these with semi-analytical results of Kirk et al. (2000). Conclusions are presented in Section 
7. 
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4 Achterberg et al. 

2 IMPORTANCE OF RELATIVISTIC SHOCKS 

In this paper we consider particle acceleration at ultra-relativistic shocks with shock velocity 
Vs = PsC and bulk Lorentz factor Fg = (1 — /3^)^^/^ ^ 1. The importance of such shocks for 
particle acceleration to energies in the EeV range (1 EeV = 10^^ eV) was pointed out by 
Hillas (1984). Quite general arguments, confirmed in our calculations below, show that the 
maximum particle energy for a particle with charge q = Ze that can be produced in a bulk 
magnetised flow on a scale Rg with velocity f3sC and magnetic field B is 

-Emax = ZeB TsPsRs ■ 

This value for -Emax is a factor Fg larger than the one obtained from the requirement that 
the particle gyration radius is roughly the size of the system, r^{E) ^ E/ZeB ~ i?s, which 
is the criterion used by Hillas (1984). This difference is due to the fact that, upstream, the 
particle typically completes only a fraction of a gyro-orbit in a regular field, corresponding 
to an angle A^ ~ 1/rs) as is explained in Section 4. This means that relativistic shocks seem 
to be the natural site of particle acceleration to extreme energies. 

First model calculations of acceleration at relativistic shocks in the context of diffusive 
shock acceleration were done by Peacock (1981). He noted that for relativistic shock speeds 
the effect of relativistic beaming on the angular distribution of accelerated particles near 
the shock becomes important. This distribution determines the energy gain and the escape 
probability per crossing cycle, the quantities that determine the slope of the distribution of 
the accelerated particles. 

This conclusion has subsequently been confirmed by semi-analytical calculations for 
shocks with Fg < 10 by Kirk & Schneider (1987a/b) and Schneider & Kirk (1989). These 
authors use an eigenfunction method to describe the particle distributions on both sides of 
the shock, and determine the slope of the spectrum from the continuity of the microscopic 
distribution function across the shock. Heavens and Drury (1988) use a similar method, and 
their results show a spectral index s ^ 2.1 — 2.3 for /3s ~ 0.98, which only depends weakly on 
the precise choice of the scattering operator employed in their calculations. The eigenfunc- 
tion method was recently extended to the case of ultra-relativistic shocks with Fg ^ 1 by 
Kirk et al. (2000). These results will serve as a independent check of the numerical results 
presented below. 

A renewed interest in ultra-relativistic shocks as sites of efficient particle acceleration 
has been sparked by key observations in two areas: 
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• The observation of Ultra-High Energy Cosmic Rays (UHECRs) in the energy range E ~ 
3-300 EeV, e.g. Bird et al. (1994); Yoshida et ah (1995) and Takeda et al. (1998). Assuming 
UHECRs are protons or hght nuclei, these particles are not confined by the magnetic field of 
our Galaxy, and should be of extragalactic origin given the fact that arrival directions do not 
cluster around the Galactic plane. Their high energy implies extraordinary circumstances in 
the production sites. 

UHECRs are observed well above the Greisen-Zatsepin-Kuz'min cut-off energy Eqzk ~ 30 
EeV, with some 10 events above 10^'^ eV. -Egzk corresponds to the energy where losses due to 
photo-pion production on the photons of the Cosmic Microwave Background become severe 
(Greisen, 1966; Zatsepin & Kuz'min, 1966). This suggests that the sources of UHECRs must 
be within a distance of ~ 50 Mpc from the Galaxy, unless UHECRs are in fact created at 
much higher energies in a 'top-down' scenario without acceleration, for instance as the decay 
products of some exotic particle (e.g. Farrar & Biermann, 1998) or some quantum-mechanical 
topological defect such as superconducting strings, e.g. Sigl et al. (1994). A review of the 
relevant observations and theoretical considerations relating to UHECRs can be found in 
Bhattacharjee & Sigl (2000). 

• The observation of the X-ray, optical and radio afterglows of Gamma Ray Bursts 
(GRBs), e.g. van Paradijs et al. (1997); Djorgovski et al. (1997); Metzger et al. (1997); 
Frail et al. (1997), which proved that GRBs originate at cosmological distances, as orig- 
inally proposed by Paczyhski (1986) and by Goodman (1986). The afterglow emission is 
believed to be Lorentz-boosted synchrotron emission from relativistic electrons (e.g. Sari, 
Piran & Narayan, 1998). These electrons must be accelerated at the external blast wave 
associated with the expanding GRB fireball, which is believed to have a bulk Lorentz factor 
Ps ^ 100 - 1000 (Cavallo & Rees, 1978; Shemi & Piran, 1990; Meszaros & Rees, 1992a/b/c, 
1993; Piran, Shemi & Narayan, 1993). The short- duration (0.01-100 seconds) gamma-ray 
fiash itself is believed to originate at internal shocks with Lorentz factor Ps ~ 2 — 10, e.g. 
Rees & Meszaros, 1994; Sari & Piran, 1995; Sari & Piran, 1997. These internal shocks are 
presumably generated by inhomogeneities in the fireball or by fluctuations in the power 
output of the 'central engine' responsible for the GRB phenomenon. A recent review of the 
GRB phenomenon and its possible implications can be found in Piran (1999). 
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6 Achterberg et al. 

A connection between these two phenomena has been suggested by Waxman (1995), who 
noted that observed UHECR flux above 100 EeV and the mean gamma-ray flux due to GRBs 
are similar, by Vietri (1995) who suggested that the production of UHECRs at the fireball 
blast wave could be extremely efficient, with an energy gain corresponding to {E{/E[) ^ 
Fg ^ 1 per crossing cycle, and by Milgrom & Usov (1995). 

Other sites where particles can be accelerated to EeV energies are the strong shocks in 
the hot spots of powerful (Fanaroff- Riley Class II) radio galaxies (Rachen & Biermann, 1993; 
Norman, Melrose & Achterberg, 1995) and in the internal shocks in the jets of Blazars. These 
shocks are expected to be mildly relativistic {Tg ^ 2 — 10). Production scenarios involving 
non-relativistic but very large shocks associated with Large Scale Structure have also been 
proposed as a source of UHECRs (Kang, Ryu & Jones, 1996). 

3 ENERGY GAIN AT RELATIVISTIC SHOCKS 

We consider a charged particle interacting with a shock with Lorentz factor Eg ^ 1 relative 
to the upstream medium. It is assumed that the upstream medium contains magnetic fields 
and possibly magnetic fluctuations which deflect or scatter charged particles in the upstream 
flow. In the case of ultra-relativistic shocks it is essential that strong magnetic fluctuations 
are present downstream so that particles can be magnetically scattered onto trajectories 
that allow them to return to the shock. This point will be considered in the discussion of 
Section 5.1. Under these circumstances a form of shock acceleration will operate with the 
same general principles as in the non-relativistic case. 

As we will show, particle acceleration near relativistic shocks is not diffusive shock ac- 
celeration because the propagation of accelerated particles near the shock, and in particular 
ahead of the shock, cannot be described as spatial diffusion. The anisotropies in the angular 
distribution of the accelerated particles is large, and the diffusion approximation for spatial 
transport does not apply. 

We consider a simple one-dimensional flow along the ^-axis. We use units where c = 1. We 
will have occasion to use three different frames of reference: the upstream rest frame (URF), 
the downstream rest frame (DRF) and the shock rest frame (SRF). Table 1 gives the notation 
used for various quantities in these frames, and the relations between these quantities as 
follow from the Lorentz transformations between these three frames of reference. A bar (^^^) 
is used to denote quantities in the DRF, and a tilde (~-) for quantities in the SRF. 
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Table 1. Fluid and particle quantities in URF, DRF and SRF for an ultra-relativistic shock 
with Ps 3> 1, and particles with 7 S> Fs, 7 2> 1- 



Quantity 




URF 






DRF 




SRF 


Shock speed 




^="^-2r| 






^s-l 







Upstream fluid speed 











^--(-I^) 




/3u = -lis 


Downstream fluid speed 




^ s 











^d~-i 


Particle energy 


E: 


= r,el (1 + AoIM) 


E 


E 


= Trol (1 - AcIM) 


E 


E 


Cosine flight direction 




_ M + Acl 

^ 1 + AelM 






(1 = 

1 - ftolM 




fJ- 


Edge of loss cone 




Sin^c = 1/Ps 






'p = 'Bs~\ 




M = 



We will assume that the shock has a speed /^s in the URF, moving in the z-direction 
towards positive z. For an ultra-relativistic shock with Fg = (1 — P^)^^'^ ^ 1 we can use 
the expansion 
1 



(3^^l- — + ^(F;^) . (3) 

If in the observer's frame the fluid moves with speed (3^ upstream, and a speed (3d down- 
stream, the velocity relevant for the shock acceleration process is the relative velocity between 
the up- and downstream gas which follows from the relativistic velocity addition law: 

This assumes that the scattering agent is passively advected by the flow on both sides of 
the shock. The velocity /9rei is the velocity of the downstream fluid seen from the URF, or 
of the upstream fluid seen from the DRF. 

For an ultra-relativistic shock with Fg ^ 1, where the magnetic fleld is dynamically 
unimportant, the jump condition relating the up- and downstream velocities in the SRF 
quickly approaches the ultra-relativistic jump condition 

^d ^ I , (5) 

regardless of the equation of state of the gas on either side of the shock. In this limit, the 
relative velocity of the up- and downstream flow satisfles (e.g. Blandford & McKee, 1976, 
see also Kirk & Duffy, 1999) 
F 1 

Trel ^-^ , Prolyl -^^ . (6) 
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8 Achterberg et al. 

This implies that in the DRF the shock moves with a velocity /?s ~ | in the positive z- 
direction. 

We will normally assume that the particles involved in the Fermi-type shock acceleration 
process are relativistic in all three reference frames, with particle Lorentz factors satisfying 
7 ^ Fs and 7 , 7^1. The particle motion relative to the shock is determined by the shock 
speed and the angle 6 between the shock normal n and the particle momentum. In the URF 
this angle is denoted by 9 and its cosine, 

fi = cosO ^ /5|| , (7) 

determines the component /3\\ of the particle's velocity /3 (in units of c) along the shock 
normal, for relativistic particles with /5 ~ 1. Particle energy and direction cosine in the 
different frames are connected by the well-known Lorentz transformations (e.g. Rybicki & 
Lightman, 1979, Ch. 4), here used in the limit p ^ E. Table 1 shows the relevant relations 
for transformations between the URF and DRF. 

These relations allow one to derive different (but equivalent) expressions for the energy 
gain which results from one shock crossing cycle. Consider a particle going through a cycle 
where it crosses the shock from upstream to downstream with energy Ei and direction cosine 
/i_d, and re-crosses into the upstream medium, after elastic scattering in the downstream 
flow where E remains constant, with an energy E{ and direction cosine /x^u- In that case 
the energy ratio as measured by an upstream observer equals 

-Ef _ 1 — /?rel/^^d /Q\ 

Ei 1 — /5rel/i_u 

One can also express the energy ratio in terms of the direction cosine fl_u measured by a 
downstream observer at the moment of re-entry into the upstream medium: 



^ = ir2 
Ei ^ ' 



Ti (1 - /?,el/i-.d) (1 + /5rel/i-.u) • (9) 

Here we have used Frei = Tg/V^. This last expression is the basis for Vietri's (1995) claim 
that particles can gain energy with Ef/Ei ~ Fg at each crossing cycle. 

In the same way, a crossing cycle where the particle leaves the downstream medium with 
direction cosine /i^u and re-enters with direction cosine /i_d after an elastic scattering or 
deflection upstream with constant E, will lead to a downstream energy ratio equal to 

Ef _ 1 + /3rel/^^u _ 1 — /3rel/i-,d /,„n 

Ei 1 + /3rel/i^d 1 — PrelfJ'-^u 

Comparing (|]) and (|lO|) shows explicitly that on average the ratio of final and initial energy 
is the same in up- and downstream rest frames. 
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3.1 Initial shock encounter 

At the first encounter with the shock, all particles with — 1 < /3/i < jSg are overtaken by 
the shock. We first consider the simplest case, that of particles almost at rest in a cold 
upstream flow so that /3 <^ 1 and Ei ^ m with m the particle rest energy. After the first 
shock crossing, the downstream energy of the particles is E = Frei m. These particles are 
immediately relativistic in the DRF if Fg ^ 1, as is assumed here. Those particles that are 
scattered back into the upstream flow have an energy equal to (see Table 1) 

Ei^lTl (l + /?,,i/i_,) m (11) 

Here we used Frei = Tg/V^. Since those particles returning upstream must have | < fi^u < 1, 
this implies that (initially) non-relativistic particles receive a boost in the upstream energy 
in the range 

r.^ < I < r.^ (12) 

at the end of their first crossing cycle. 

If the upstream medium contains a pre-existing population of relativistic particles, dis- 
tributed isotropically in the URF, one must average Eqn. (|^) over all direction cosines 
—1 < yU_d < j3s with a weight proportional to the particle flux oc /5s — /i_d across the shock 
into the downstream medium. This yields an average initial boost equal to 

F^ (l+/?,el/^^,) , (13) 

which lies in the range 

9 s ^ ^ — 3-^ s ■ K^^J 

Here we assume that /i^u is uncorrelated with /x^d- 

We conclude that in all cases the energy gain at the first shock encounter is of order 
Ef/Ei ~ Fg for those particles which return to the upstream medium, so that Vietri's (1995) 
suggestion applies at the first shock encounter. This provides a natural injection process 
for further particle acceleration by repeated shock crossings at ultra-relativistic shocks. We 
return to this point in Section 6.3. 

This conclusion can also be reached by considering the Rankine-Hugoniot relations for 
an ultra-relativistic shock (e.g. Blandford & McKee, 1976) since, from a microscopic point of 
view, the situation is similar: in a shock the kinetic energy of the upstream flow must be 'ran- 
domized' so that a hot downstream state arises behind the shock. The proper downstream 
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10 Achterberg et al. 

energy density e, pressure P and number density n of the downstream flow are related to 
the corresponding upstream quantities by 

e^W^ 2Tl (e + P) , n ^ 2^2 Tg n . 

The typical 'thermal' energy per particle in the downstream is 

n \ n / 

In the case of a cold upstream medium (P <^ e ~ nm) one finds E^\^ k. T^eym, with T-^^x ~ 
Fs/v^ for an ultra-relativistic shock, while in the case of a relativistically hot upstream fluid 
(P = e/3) one finds ii^th = | Trei E, with E = e/n the mean upstream energy per particle. A 
Lorentz-transformation back to the URF adds another factor ~ Fg, so an upstream observer 
assigns these particles in both cases an energy Ef ^ F^ E^ with Ei the initial upstream 
energy. 

3.2 Kinematical constraints on the energy gain at relativistic shocks 

We now consider true Fermi-type acceleration at an ultra-relativistic shock of particles that 
have been downstream at least once. Particles re-crossing the shock into the upstream region 
satisfy at the moment of shock crossing 

7I>;^^| , yU>/3s (15) 

for the down- and upstream particle flight direction respectively. These relations define the 
downstream loss cone with a corresponding opening angle 6c around the shock normal in 

the upstream frame. This opening angle follows from 

1 



sin9c = sJl-(3!- . 

J- s 

As 6c is small for Fg ^ 1 we can use sin 6c ~ 6c to write the condition ([15|) for crossing the 
shock into the upstream flow in terms of upstream variables as 

^^u < ^c ~ ^ . (16) 

J- s 

Particles that have just entered the upstream region reside within this loss cone. Upstream 
deflection or scattering must change the upstream flight angle to a value 6 > 6c before a 
new shock crossing cycle can begin. 

As we will show below, all plausible scattering and deflection mechanisms are only able to 
change the angle 6 by an amount \A6\ ~ 6c before particles are again overtaken by the shock. 
As a result, the upstream angular distribution of these particles is confined to a cone with 

© 0000 RAS, MNRAS 000, 000-000 



Particle acceleration 11 

opening angle A9 ~ 2/rs ^ 1 around the shock norniah an effect of relativistic beaming. 
For ultra-relativistic shocks the small angle approximation, sin 6' k, 0, is an excellent one for 
all particles in the upstream region which have interacted with the shock at least once, in 
particular those about to be overtaken by the shock with 1 < rs6'_d <2. 

Consider a particle performing a crossing cycle from upstream to downstream and back. 
The angles in this crossing cycle satisfy 

^^u < ^ < ^^d . (17) 

If we use the small angle approximation /z = 1 — |^^ in expression (|^) for the energy ratio 
in the URF at the end of the cycle, together with expansion (^ for /3rei, one finds (Gallant 
& Achterberg, 1999): 

Ei^ 2 + T^,9l^ ■ ^ ^ 

Because of the inequality (|1^) this ratio is always larger than unity. But on average (as we 
will demonstrate explicitly in our simulations) its value is never much larger than Ef/Ei ~ 2 
since both the angular factors in this expression remain of order unity. A similar conclusion 
holds for the energy ratio Ei/E\ seen by a downstream observer in a cycle where the particle 
moves from downstream to upstream and back, as is readily seen from Eqn. (p!0|). 

This shows that Vietri's (1995) suggestion does not apply to these particles. If one uses 
expression (H) for the energy ratio in a cycle starting in the upstream region, the factor 

l-Ael/i-.d^^ + ^ (19) 

is always of order F"^ due to relativistic beaming, while the factor | < (1 + /5rei/i_u) < 2 is 
always of order unity. Once again one concludes that the value of Ef/Ei niust be of order 
unity. 

4 DEFLECTION AND SCATTERING 

As stated above, the process of diffusive shock acceleration in the non-relativistic case relies 
on scattering on magnetic irregularities to confine the accelerated particles near the shock 
for some time. For the case of ultra-relativistic shocks we will consider two cases in some 
detail: the case of regular deflection upstream and strong scattering downstream, and the 
case where strong scattering operates on both sides of the shock. 
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4.1 Regular deflection 

The simplest way upstream particles will leave the loss cone is through deflection by the 
Lorentz force in a uniform upstream magnetic field. Without loss of generality we will choose 
the shock normal in the z-direction, while the magnetic field is assumed to lie in the x — z 
plane: B = [B^ , Q , B\\). 

The equation of motion for an ultra-relativistic particle of charge q = Ze and energy E 
with gyration frequency VL^ = ZeBc/E in a magnetic field B = B b reads 

^ = fig(/3x6). (20) 

This equation can be solved approximately using the fact that 

A = i-o(r:^) , I/3J ^ ^ = o(r;i) « 1 

at all times when the particle is upstream. Here (3\\ ~ /x is the component of the particle 
velocity along the shock normal, and (5^ = {Px , Py ,0) the velocity component in the plane 
of the shock, all measured in the upstream frame. Unless the field is almost along the shock 
normal so that 

one can solve (^) by iteration, with only Q± = ZeB±c/E appearing in the equations to 
leading order in e. We only consider the case e -C 1, covering almost all possible field 
orientations except for a fraction ~ l/rg of the solid angle. We assume that the particles 
are ultra-relativistic so that we can put f3 = Jf^u + (3j_ = 1, and choose ZeB± > without 
loss of generality. 

Under these assumptions one finds for a particle that crosses the shock at t = and 
z = with initial momentum direction corresponding to /3(0) = {Pxi , Pyi , PzC)'- 



PAt) = P: 



XI } 



Py{t) ^ Pyi + n^t, (22) 



z{t) ^ p.it-lPyi{n^tf~l {n^tf . 



The expression for z{t) = J^dt' Pz{t') follows from /J^ ^ 1 — \P']_- We have consistently 
neglected all terms of order l/Fs and e with respect to unity in the expressions for P^ 
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and (3y. The time t^ of the next encounter with the shock, the upstream residence time, is 
obtained by putting z{tu) = z^ = Pstu, which yields: 



^±iu — 2 



9/^S + 12 (^ - /5io 



ni/2 

— 2h'yi 



i/3,i . (23) 



The corresponding values of Pyit^) = Pyt and the position Zs can be found from substituting 
this expression in Eqn. (|22|) . In particular we find Pxi = Pxi and 



/^.f = -^ + \/fI-3/5ii-4/?i, (24) 

c.f. Gallant & Achterberg (1999). The angle with respect to the shock normal follows from 



(31 + Pi = sin^ e^e'. (25) 

In figure 1 we give a graphical representation of the position of the particle in the (3x-Py 
plane when it is overtaken by the shock at t = tu, as a function of the initial conditions {(3xi, 
(3yi) at the moment the particle enters the upstream region. 

The largest defiections for given (3yi occur when the particles enter the upstream region 
with (5x\ = 0. In that case we have in the small angle approximation 



and particle trajectories remain in the y — z plane. In figure 2 we give the value of (3yf = ^_d, 
the energy gain and the upstream residence time as rgfi^tu- The largest value for 6^^ occurs 
when Fg (3yi = —1: then Fg 6^^ = 2. However, the largest energy change from Eqn. (10), 
corresponding to Ef/Ei ^ 2.62, occurs for TsPyi ~ —0.27. 

These results set the following limits on the angle ^^d and its downstream equivalent 
^^d = cos~^(/i^d) with respect to the shock normal when particles enter the downstream 
medium after regular defiection by an upstream magnetic field: 

l<F3e^d<2 ^^ |>/i^d>-|. (26) 

4.2 Scattering 

We now consider the circumstances under which scattering rather than defiection by an 
ambient upstream field changes the upstream particle momentum. In the case of ultra- 
relativistic shocks the small opening angle 6^ ~ l/Fg of the loss cone strongly limits the 
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initial and final orientation 



- r„^i, = 0, 0.1 0.9, 0.95 



■d 
Id 




Figure 1. The location of a particle in the Vsf^x-^sPy plane when it is overtaken by the shock as a function of its location when 
it enters the upstream region. Shown is the location of particles entering on concentric circles within the loss cone centered 
on the shock normal. The dashed circle with radius Fg/Sx = 1 corresponds to the edge of the loss cone. The concentric circles 
have radii corresponding to Fs/Sxi equal to 0, 0.1, 0.2, ■ ■ ■, 0.9 and 0.95 (the dot and the thick solid circles). These circles map 
onto the kidney-shaped curves along lines of constant fix- The origin /3^i = (3y\ = 0, corresponding to a particle entering the 
upstream flow along the shock normal, maps to /3j.f = 0, TsPyi = a/S. This is indicated by the thick dots. The larger changes 
in the angle 9 ft: |/3x| occurs for particles with Py\ < and fi-^i = 0, with the largest change ocurring for V^Pyi = —1- These 
particles must turn through most of the loss cone before they leave it and are overtaken by the shock. In contrast, particles 
with /3yi > change their orientation by a relatively small amount, leaving the loss cone almost immediately. This figure is for 
ZeB^ > 0. If the perpendicular field has the opposite sign, this figure should be mirrored with respect to the /3i-axis. 

regime in which the direction of the upstream particle momentum can exhibit diffusive 
behaviour. 

As a simple model for scattering of a charged particle by waves or random magnetic fields 
consider a process where the angle between shock normal and particle momentum changes 
at a constant rate, whose sign is chosen at random, during a correlation time tc ■, 



9{t + At) = 9{t) ± 9,^, ( —] for At < t, 



(27) 



For t ^ tc this process corresponds to a random walk in 9 with step size 6'rms and corre- 
sponding diffusion coefficient 
9' 



Va 



?2 
rms 



An upstream particle will only enter the regime where the momentum direction diffuses if 
it is not swept up by the shock within a correlation time tc- Since the regular motion within 
a correlation time as decribed by (|27| ) is formally equivalent with the case considered in 
the previous Section with /^j, = 0, |/?y| = 6* and VLi_ — > 6'rms /^c, the capture time t^ in the 
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final angle, turning angle and energy ratio 
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Figure 2. The final angle 8^^ ~ /3yf, the upstream residence time, plotted as FsHxtu, and the downstream energy gain 
E[/Ei = (i?f/Ei)dwn in the downstream-upstream-downstream cycle starting with S_u = \f3yi\ and ending with 6_d = /3y(, all 
as a function of f3yi with respect to the shock normal at the moment the particle enters the upstream fiow. Curves are shown 
for those particles moving in the y-z plane, i.e. f3^i = f3^{ = 0. 



regular regime for a particle crossing the shock with an angle < 0_u < V^s follows from 
(p3|) as 

3 



9- 



^-^euT^^^u<^ 



(29) 



where the minus (plus) sign applies in the half space Pyi = 9^^ and Pyi = —9^^ respectively. 
The maximum value of t^ occurs for TsPyi = — 1 (see figure 2). This equation has a solution 
with tu < tc for all allowed values of 9^^ provided 

^rms >^ = S9,. (30) 

J- s 

If this condition on the stepsize of the random walk in 9 is satisfied, upstream particles are 
overtaken by the shock before the direction of their momentum starts to diffuse. In that 
case, all the results for regular deflection apply. 

If condition ( PPf ) is not satisfied, diffusion will start before the particle is overtaken. The 
average position at time t of particles entering the upstream medium at t = at an angle 
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0^^ with respect to the shock normal can be calculated as (see Appendix B of Achterberg 
et al. (1999) using units with c = 1) 

i^z{t)) = {l-\dl,)t-Vet\ (31) 

assuming that Dgt <^ 1 but t ^ tc- The average value of 6'^ satisfies in the same limit 

e\t)) = 9i^^ + 4:Vgt. (32) 

The diffusive change in 6"^ differs from the naive estimate, {0'^{t)) — ^^^ = ^Dgt, due to the 
effect of dynamical friction. The typical time t^ at which a particle is overtaken by the shock 
can be estimated by putting ^s = Pstu = {z(tu))- This yields the relation 

^otu ~ ^ (i - r^'^u) , (33) 

and allows us to estimate the angle at which the particle re-crosses into the downstream 
medium as (Gallant & Achterberg, 1999) 

Oi,-(0'iQ) = l,-Oi^. (34) 



Y2 

s 

These relations are only approximate since they neglect the correlation between the shock 
re-crossing time and the re-crossing angle. Note that they do satisfy the intuitive requirement 
that in this regime particles entering the upstream region near the edge of the loss cone, 
6'^u~ l/fs, are recaptured almost immediately by the shock just outside the loss cone with 

^^d>i/rs. 

The typical downstream ratio of the final and initial energy in a crossing cycle where the 
particle moves from downstream to upstream and back follows from Eqns. ([I0|) and (0) as 

This downstream energy ratio falls steadily from Ef/Ei = 2 at 9^^ = to Ef/E^ = 1 at 
the edge of the loss cone where TgO^u = 1- Note however that this is an estimate based on 
averages, so that invidual particles can have a larger energy gain. 

4.3 Scattering agents: waves and random magnetic fields 

In diffusive acceleration near non-relativistic shocks one usually relies on gyro-resonant scat- 
tering of the particles by low-frequency MHD waves (Alfven waves) supported by a nearly 
uniform magnetic field (e.g. Wentzel, 1974). The magnetic field perturbation 6B in these 
waves satisfies 6B ■ B = 0. In this process particles interact mainly with those waves that 
satisfy the resonance condition 
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uj{k) — k^\V\\ ± fig = , 

where uj{k) is the wave frequency and k is the wave vector, with fcy the component along 
the mean magnetic field. For low-frequency waves this usually means that their wavelength 
must be of order of the particle gyration radius. 

In the rather extreme circumstances considered here, this process will only be of impor- 
tance if B^ <^ BjTg so that defiection by the large-scale field can be neglected to lowest 
order. Note that this is only the case for a small fraction ^ l/^s of possible field orien- 
tations. We will consider the limiting case B^ = and -By = B. The resonant scattering 
process itself corresponds to a random walk in 6 with 

6*^8 ~ ( ~B" ) ' ^c ~ 27r/fig , 

with 6B the typical amplitude of the waves at resonant wavelengths. If (|30|) is satisfied, 
corresponding to a wave amplitude 

5B > 3B/T, , (36) 

upstream particles will never enter into the diffusive regime. In that case they spend less 
than a gyroperiod upstream, and the resonance between waves and particle can not be 
achieved. In this regime the non-resonant 'sloshing motion' of the particles in the magnetic 
field of waves with wavelength A>rg, which has an amplitude A6 ^ 6B/B, is sufficient to 
defiect them out of the loss cone so that they are overtaken. This means that only for very 
small wave amplitudes does true diffusion take place upstream, with a diffusion coefficient 
(Wentzel, 1974) 

Ve ~ fig (^) . (37) 

If no large-scale upstream magnetic field is present, similar conclusions hold. Consider 
a random field with rms amplitude Bj- and coherence length ic- The simplest realization of 
this situation is a collection of cells of size 2ic with a random orientation of the magnetic 
field B and a field strength \B\ = B^.. The typical gyration radius of a particle in such a 
field is 

r^E) = ^ . (38) 

Charged particles traversing such a random field exhibit a random walk in the momentum 
direction of the sort discussed in the previous Section, provided ic <^ rg{E). This random 
walk can be characterised by (Achterberg et al. 1999): 
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4 24 

6'rms ~ o — , ic ~ • (39) 

■^ Tg C 

The deflection of a particle in a single magnetic cell will be sufficient to lead to a shock 
crossing if (^) is satisfied, which in this case reduces to 

4 > ! ^ . (40) 

J- s 

This condition is equivalent to 

E <E, = I ZeB, T, 4 . (41) 

Only particles with E > E^ will be in the diffusive regime with diffusion coefficient 



cL 



V, = ^_. (42) 

These calculations show that upstream, true diffusion of the particle flight direction occurs 
only if the large-scale upstream magnetic fleld (if present) is closely aligned with the shock 
normal, and if the deflection by waves or magnetic cells occurs in sufficiently small steps. 
If this is not the case, particle sloshing motion in the waves or deflection within a single 
magnetic cell is already sufficient to deflect a particle out of the loss cone so that a new 
shock crossing results. 

In the downstream medium, on the other hand, particles have to turn through a large 
angle of order A^ ~ 1 in order to re-cross the shock, and strong diffusion in the angle 9 is 
needed. 



5 CYCLE TIME AND MAXIMUM ENERGY 

The maximum energy that can be achieved in Fermi-type acceleration at an ultra-relativistic 
shock can be calculated in the same way as in the non-relativistic case: acceleration proceeds 
until the energy gain per cycle becomes equal to the energy losses incurred over a shock 
crossing cycle. Since the typical energy gain per cycle at an ultra-relativistic shock satisfies 

i\E = E{ — E\ ~ i?i , 

the acceleration time scale and the cycle time, which consists of the time a particle resides 
in the up- and downstream region, tcy{E) = t^ + t^, are essentially the same. Since the 
relative speed between the URF and DRF is relativistic, one has to be careful about the 
frame used in calculating the two contributions to tcy In particular, if a particle of energy E 
in the DRF spends a time td in that frame before re-crossing the shock, the corresponding 
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upstream energy is -E ~ Fg E, and the downstream residence time, separating two events at 
the shock, is equal to tj ~ Tg td- This means that the cycle time measured by an upstream 
observer can be written as 

t,y{E)r^t^{E) + T,UE/T,). (43) 

5.1 Cycle time for regular deflection 

Let us first consider the upstream residence time in the case of defiection by a large-scale 
magnetic field. In that case the particle momentum must typically turn through an angle 
A^ ~ 1/Ts before the particle is overtaken once again by the shock. This implies 

tn[E) = „ „ ^ , (44) 

in accordance with the results of Section 4.1. 

In the downstream medium we must appeal to some form of scattering. If the downstream 
field consists simply of the compressed large-scale upstream field resulting from the shock 
jump conditions, where the perpendicular field component is amplified according to B± ^ 
2^/2TsB± while the parallel field component remains the same, the downstream magnetic 
field will be almost completely aligned with the shock surface. In that case shock acceleration 
can only result if there is efficient cross-field diffusion of particles ( Jokipii, 1987; Achterberg 
& Ball, 1994) so that particles can catch up with the shock, which moves with a speed ~ c/3 
with respect to the downstream medium. This diffusion is presumably due to (shock-induced) 
wave turbulence. 

This requirement explains why the simulations of Ballard and Heavens (1992), who follow 
exact particle orbits in a random magnetic field, find much steeper spectra for particles 
interacting with a relativistic shock {f3s ~ 0.98). Without a stochastic process which allows 
for a rapid decorrelation between particle and magnetic field in the downstream region, only 
particles which are located on a field line with a favorable geometry can re-cross the shock. 
This reduces the average return probability which leads to a steepening of the spectrum. 

As long as the turbulence is not so strong that the downstream field geometry becomes 
completely chaotic, such a situation can only be realised if the cross-field diffusion coeffi- 
cient D± is close to the maximum possible value: Bohm diffusion with diffusion coefficient 
-Db = crj/3, where particles randomly move across the field with a step size of order the 
particle gyration radius every gyration period. This follows from the requirement that the 
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diffusive stepsize A'z along the shock normal, which downstream is almost perpendicular to 
the downstream field, must satisfy (Achterberg & Ball, 1994) 

/S.Z > cP^T ^ cr/3 , 

with r the time interval between diffusive steps. For cross-field diffusion in the weak turbu- 
lence limit one has Az ~ f^, with f^ the downstream gyroradius. The above condition then 
corresponds to r < 3r^/c and to a diffusion coefficient which must satisfy 

-2 



' 'b 



2r 6 "^ 



This is the case on which we will base our estimate below (Eqn. ^6D for the downstream 
residence time. 

If the downstream field is completely chaotic, leading to quasi-isotropic diffusion of the 
momentum direction downstream, particles in general will find their way back to the shock, 
as our simulations show. We will show results for the typical downstream residence time for 
this case in Section 6. 

The typical residence time of a particle in the downstream medium is (Drury, 1983) 
_ _^ AD^ AE [DA 

'^<^) - wr ^eSi^Av;) ■ («' 

Simple shock compression leads to a downstream magnetic field with B± = 2\/2TsB±, so 
scaling law (^) gives the corresponding time in the URF for Bohm diffusion {D± ^ Db) as 

Ms) = r.l.(B/r.) = ^^. (46) 



Comparing this with expression (44) one sees that td ~ tu at a given energy. If the cross- 
field diffusion is strongly enhanced by the effect of field-line wandering (Achterberg & Ball, 
1994) in a turbulent downstream flow, the downstream residence time could be increased 
significantly so that the particle spends most of its time downstream. On the other hand, 
if turbulence is strongly enhanced downstream so that scattering proceeds at a rate with 
effective collision frequency u^ = l/t^ ^ fig, particles are no longer effectively tied to field 
lines. In that case particles will diffuse almost isotropically downstream, and will spend most 
of a shock crossing cycle in the upstream flow. 

5.2 Cycle time for upstream momentum diffusion 

If there is no large-scale field, but rather a randomly oriented field with coherence length 
ic, and ii E> ZeB^Tsic so that the momentum direction diffuses between shock crossings, 
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the upstream residence time is roughly the time it takes the momentum direction to diffuse 



through an angle ~ l/Fg in the random field. Using (^^ one finds 



-^-^-rk-i!k.i:^y 



s^c 



The extra factor with respect to the regular case is always larger than unity because of 
condition £c ~ ^g/Fg for the diffusion regime to apply upstream. 

A similar calculation downstream, where particles must typically turn through an angle 
A^ ~ 1, yields: 

The amount of additional field amplification downstream above shock compression can be 
parameterised by the ratio 

Using 

-,{E)^r^{E)/{2V2Tli^), 

one can write the corresponding time measured by an upstream observer as 

tJE) = TJAE/T,) ^ — lT^\ . (50) 

Downstream the particles will diffuse provided 4 < r^- If downstream field amplification 
is large, ^b ^ 1, as is often assumed in afterglow models for Gamma Ray Bursts, the 
downstream residence time would typically be less than the upstream residence time. 

5.3 Maximum energy estimates 

The above estimates for the cycle time allow one to calculate the maximum energy that can 
be achieved in Fermi-type acceleration at an ultra-relativistic shock, generalizing the results 
of Lagage and Cesarsky (1983), Drury (1983) and Heavens (1984) to the relativistic case. 
Consider a propagating spherical shock with radius Rs- The finite age of the shock, tg ^ Rs/c, 
limits the energy of shock-accelerated particles through the requirement tcy{E) < Rs/c. The 
best possible case arises when td <^ tu- For regular deflection upstream one then finds 

E^ax ~ ZeB^ F,i?, . (51) 

If the particle momentum starts to diffuse in a random upstream field, which requires 
E>E^r^ ZeB^ Fs^c, the maximum energy equals 
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E^^^^ZeBrV, (RJjsf^ . (52) 

This expression for -Emax is only relevant if Rs > ic- 

Similar limits on the maximum energy apply to a standing shock in a spherical wind or 
in a jet with constant opening angle. Expansion losses of relativistic particles in a relativistic 
flow satisfy (e.g. Webb, 1985) 

(53) 



'dE'\ E' 



dt' / 

/ exp 



dV 



Here E' and t' are the particle energy and time measured in the fluid rest frame, and V{x, t) 
and r(a;, t) are the fluid velocity and associated Lorentz factor in the observer's frame. This 
result applies if frequent scattering keeps the particle distribution nearly isotropic in the 
fluid rest frame. In that case one can use the Lorentz-invariance of the energy loss rate, 
dE'/dt' = dE/dt (e.g. Rybicki & Lightman, 1979, Ch. 4), together with E ~ TE'. These 
relations then yield the expansion losses of relativistic particles {E ^ pc) in the observer's 
frame: 



dE\ 


E 


dr 


dt 


3r 


dt 



exp 



V ■ {TV) 



Therefore, expansion losses in a flow with typical velocity Vs and size Rs are of order 
fdE\ /K 



(54) 



dt / Vi?s ' 

\ / exp '^ 

and involve the same dynamical timescale ts ~ Rs/Vs- Particles will no longer have a net 
energy gain in a crossing cycle if tcy{E) > ts- 

The above estimates for -Emax neglect radiation losses. For electrons (or positrons) these 
losses can limit the maximum energy to lower values. This was discussed for non-relativistic 
shocks by Blandford(1977) and Heavens (1984). The energy loss rate for an electron with 
E ^ rricC^ scales as 

f = -|axcf/. f-^)'. (55) 

dt •^ KmeC^J 

Here o"t is the Thomson cross section and 
B' 

2 



Ut = f/rad + 9 - — • 



TT 

The first term in Ut is the energy density of the radiation field (assumed to be isotropic) 
leading to inverse Compton losses. The second term corresponds to synchrotron losses in 
the ambient magnetic field. Since all upstream particles are closely aligned with the shock 
normal only the field component Bj^ contributes to the synchrotron losses. 
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Consider a particle entering the upstream medium at t = with an energy Ei. Its energy 
decays due to radiation losses as 

m = — f— , (56) 

with 

_ SmgC^ f mcC^\ 

the loss time at energy E-^. Assuming regular deflection upstream, the angle 6 between 

particle momentum and shock normal changes according to 

d^ ZeB^c ^ / t\ .^. 

with fii = ZeB±c/Ei. The time to turn through an angle A6 is 



t Tgy 



1 2Ag ^ 



and the electron energy is reduced in that time to 
E, 



E 



2A9 
1 



(58) 



\ ' 'i ''^sy 

Since shock crossings typically double the particle energy in each cycle, and since Ad ~ l/Fg, 
upstream losses will lead to a net energy loss in a cycle if 

Q.-r >i/r 

which corresponds to 



E > Esu = rrieC 



,2 



'3e5,r.\i/2 



Losses in the downstream leg of the crossing cycle will become important if the downstream 

residence time and the synchrotron loss time are roughly equal. Assuming Bohm diffusion 

(Eqn. |45| ) this condition reads 

4:E 3meC^ (meC^\ 

ZeB±c ~ AaTcUt \ E ) ' 
which corresponds to an energy 

- _ 2 / 3e5| 

E^ E,d = m^c' I - 



1/2 



\16crTf/t/ 

in the DRF. In the URF these particles have an energy 

E,d ~ Ts E,A = Ts m,c^ (-^^] ■ (60) 
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If Compton losses can be neglected (t/t ~ 3i?^/167r) one has 



2 f^T^eT,\'/^ r. ^.r.n ^s 



Esu ~ m,c' —-^ , ^sd ~ 0.3 _ . 

Here we have used definition (|49|) in the second estimate, where the numerical constant 
is l/(8v^)^''^- This shows that for ,^b ~ 1, up- and downstream losses are roughly equally 
important, but that the downstream losses will determine the maximum energy if significant 
field amplification occurs in the downstream fiow. If losses are important, the maximum 
energy of shock-accelerated particles will be 

-E'max = niin {Eg^ , Esd) ■ 

6 SIMULATIONS 

We have performed a number of numerical simulations of the acceleration process using the 
following assumptions: 

• The defiection mechanism in the upstream region is either pure scattering, or pure 
defiection by a homogeneous upstream field B±] 

• In the downstream region there is strong (diffusive) scattering and no regular deflection; 

• Up- and downstream fluid states are connected by the ultra-relativistic jump condition 

PuPd 31 

• Radiation losses are neglected. 

Our simulations employ the method of Stochastic Differential Equations introduced by Ito 
(e.g. Gardiner, 1983; Oksendal, 1991; Achterberg & Kriills, 1992). We follow particles in the 
shock rest frame, advancing the particle position z along the shock normal at each time step 

as 

Az = fu,d (n^ - /5u,d) Asu,d • (61) 

Here n^ is the component of the unit vector h = p/\p\ along the shock normal, as measured in 
the URF or DRF respectively, and Asu,d = c Atu^ the path length increase of the particle in 
the upstream or downstream fluid rest frame, which is taken to be a constant. The quantities 
ru,d and /5u,d are the up- and downstream fluid Lorentz factor and fluid speed in units of c 
in the SRF. This expression follows straightforwardly from the Lorentz transformation for 
position between the SRF and the URF or DRF. The shock is located at ^ = 0, and when 
a shock crossing is detected the particle momentum and energy is transformed between the 
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URF and DRF using exact Lorentz transformations. Particles are then propagated further, 
starting at the shock. 

To describe scattering, the particle flight direction n in the URF and DRF is advanced 
according to (e.g. Achterberg et al., 1999) 



n{s + As) = V 1 - |Anst|^ fi{s) + An^t . (62) 

In this expression Arigt is the stochastic change in the orientation of the flight direction which 
results from scattering by random fields or waves. Expressed in terms of the corresponding 
angular diffusion coefficient T>q, this change in the flight direction is given by 



An,t = pVe (As/c) (^lei + ^62) , (63) 

where ^1 and ^2 are two independent unit Wiener processes, drawn at each step from a 
Gaussian distribution with unit dispersion so that they obey the simple rules (^1) = (^2) = 
O5 {ii) = {&) = 1 a^d (^1^2) = 0, where the brackets indicate an average over many 
independent steps. The two (arbitrary but mutually orthogonal) unit vectors 61^2 are chosen 
in the plane perpendicular to n so that Aust • n = 0. Note that the norm n h = 1 is 
preserved identically. It can be shown that this gives an excellent approximation for isotropic 
scattering (i.e. T>g independent oi 6) for |Anst| ^ 1. 

Upstream, the value of As is chosen in such a way that the typical diffusive step A^st 
in the flight direction satisfies FgA^st < 0.1. The downstream value As is chosen in such a 
way that it matches the upstream value after the appropriate Lorentz transformation, until 
the particle has diffused a distance of a few mean free paths behind the shock. Then As is 
typically chosen ten times larger. 

In the case of regular deflection by a homogeneous field in the upstream medium the 
analytical solution (|^) is used to calculate the orientation angle 6^d of the particle momen- 
tum vector at the moment the particle recrosses the shock, given the flight direction 6^-^ 
with which the particle enters the upstream region. 

By running many, statistically independent, realisations of this prescription and recording 
the particle energy and momentum direction at each shock encounter, one can construct the 
particle distribution in momentum space at the shock which results from repeated shock 
encounters. This corresponds to the situation where a steady-state is reached in the shock 
frame. 

We also use a particle splitting method, where the loss of particles from the acceleration 
process due to the large value of Pesc is compensated by the introduction of additional 
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Table 2. The spectral slope s, the mean ratio of initial and final energy < E(/Ei > in a crossing cycle upstream > downstream 

> upstream, and the return probability Prot for the case of strong scattering on both sides of the shock. 



Fa S < Ei/Ei > V,ot 

10 2.230 ± 0.012 1.97 ± 0.01 0.435 ± 0.004 

10^ 2.219 ± 0.004 1.97 ± 0.01 0.439 ± 0.006 

103 2.213 ± 0.003 1.93 ± 0.02 0.437 ± 0.006 



particles, so that the effects of Poisson noise in the simulations is minimized. Typically 
particles are split every time they complete five crossing cycles, by keeping track of the 
number ridu of downstream to upstream crossings. In this splitting method, a particle is 
split in into A^ copy particles, each with weight W = 1/N. We typically use A^ = 25 — 100. 
These copy particles start their evolution with the same initial conditions (nju, momentum 
and position) as the parent particle. Thereafter, all copy particles evolve independenty, 
according to the stochastic algorithm outlined above. 

The splitting procedure has no noticable effect on the particle distribution. This is most 
easily seen by looking at the distribution of shock crossings, as this quantity determines the 
particle splitting. The distribution is a featureless power law, with the slope set by the escape 
probability Pesc (see figure 3). Neither the position of the splitting boundaries in terms of the 
number of shock crossings, nor the number of particles created at each boundary influences 
this result in any noticable way. The same holds for the energy distribution (figure 4) which, 
apart from the transient effect of the injection conditions at the beginning of the simulation, 
is once again a featureless power law in energy. 

An absorbing boundary is placed sufficiently far downstream (at a distance from the 
shock at least 10 x the scattering mean free path) to remove particles from the simulation 
which have a vanishingly small chance of ever returning to the shock. 

6.1 Case of pure scattering 

We first consider the case where particles are scattered on both sides of the shock. For this 
case the semi-analytical results from Kirk et al. (2000) can serve as a basis for comparison. 
Table 2 gives the results for a number of simulations at different values of Fg. 

The results presented in Sections 3 and 4 of this paper are the leading terms in a formal 
expansion in powers of l/F^. This means that all the simulations for Fg = 10, 100 and 1000 
should give identical results for the slope of the distribution, the mean energy gain per cycle 
and the return probability within about 1%. Table 2 shows that this is indeed the case. 
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Figure 3. The distribution of the number of crossings n^u from the downstream region into the upstream region. The distribu- 
tion is a featureless power-law, as should be the case for particles with an energy-independent return probabihty Pret = 1 — 'Pesc- 
The distribution shows no sign of an influence of the particle splitting employed to diminish the Poisson noise, which uses n^u 
in the sphtting criterion. 

We can compare the slope s of the energy distribution of the accelerated particles as found 
in these simulations with the analytical results of Kirk et al. (2000), and the simulations 
of Bednarz & Ostrowski (1998). In our simulations s is obtained from a least square fitting 
method on the distribution of accelerated particles, at energies exceeding at least 30 x the 
injection energy so that there is no 'memory' of the initial conditions at injection. The 
various values agree within quoted errors: 

2.22 ± 0.01 these simulations; 

2.2 Bednarz & Ostrowski (1998); 

2.23 ±0.01 Kirket al. (2000); 

As a consistency check on the value of s obtained from the simulations, we can also use 
relation (|l|) to calculate it from the return probability Pj-et and from the mean energy gain per 
crossing {Ef/Ei) obtained in these simulations. These are obtained by recording the number 
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Figure 4. The distribution of particles as a function of the upstream energy gain (Ef/E{)^, for particles interacting with a 
shock with Fg = 1000. The particles are injected upstream at energy E[ with their flight direction distributed randomly within 
the loss cone 9 < l/Fg. A featureless power-law establishes itself after a few crossings, signalling that the memory of the initial 
conditions has been erased. 

of crossing cycles a particle completes, and the energy gain per crossing. Substituting the 
values from Table 2 we find 

ln(l/Pret) 



1 + 



2.23 ±0.02 



In {Ef/E,) 

Figures 5 and 6 shows a more detailed comparison between the simulation results and the 
semi-analytical results of Kirk et al. (2000). 

Figure 5 shows the angular distribution of particles at the shock as a function of the cosine 
of the angle between particle momentum and shock normal, and Figure 6 shows the flux 
distribution. At the position of the edge of the loss cone (cos(6'd) = |) there is an unphysical 
spike in the angular distribution. This spike is a discretization error, resulting from the fact 
that particles take discrete steps in angle, c.f. Eqn. (|63D. This allows particles to be 'trapped' 
near the edge of the loss cone since they can re-cross the shock after only one step. In that 
limit, scattering is not well-represented by our method. This angular distribution is in fact 
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Figure 5. Tiie angular distribution of particles at the shock in the downstream rest frame. The angle 8^ is the angle between 
the particle momentum and the shock normal, with ^^ = corresponding to a particle crossing the shock along this normal into 
the upstream medium. The distribution is normalized to a unit integral. The histogram gives the results from our simulations, 
and the smooth curve shows the semi-analytical result from the eigenfunction approach. The downstream loss cone corresponds 
to cos 8^ > /3g ai TT. This curve is for particles with an energy E > 30 X -Einji with Ei^j the energy of a particle when it starts 
true (Fermi-type) shock acceleration. This ensures that the distribution has relaxed so that the initial conditions no longer 
influence the shape of the distribution. 

reconstructed from the particle flux across the shock, which is the quantity recorded in our 
simulations. 

We feel confident that this peak does not influence the results for the slope s of the 
spectrum of accelerated particles. The relevant parameters for the acceleration of particles, 
such as the mean energy ratio {Ef/E[) and the return probability Pret? are in fact flux 
averages. Particles near the edge of the loss cone contribute a negligible flux, and the spike 
is absent from the flux distribution of Figure 6. The flux distribution obtained from the 
simulations agrees well with the analytical results using the eigenfunction method described 
in Kirk et al. (2000). 

Figure 7 shows the distribution of the upstream residence time ty,{E) and the downstream 
residence time td for Fg = 1000, both measured in the respective fluid rest frames for particles 
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Figure 6. The distribution of the absolute value of the flux J^ across the shock as a function of cos(Sd) in downstream rest 
frame variables. The histogram gives the results from our simulations, and the smooth curve is the semi-analytical result from 
the eigenfunction approach of Kirk et al. (2000). The flux vanishes at the edge of the downstream loss cone, cos 9^ = /3s ~ ^i 
where particles move exactly along the face of the shock. This curve is for particles with an energy B > 30 X E\^y where the 
initial conditions no longer influence the distribution. 

with i? > 30 X i^inj. The distribution is plotted as tu (dA^/dtu) and its downstream equivalent 
as a function of the dimensionless variable 
FgDetu upstream 

r = < 

"Detd downstream 

According to our discussion in Section 5 both distributions should peak near the average 
value of the residence time which should lie in the range r = 0{1). 



Both distributions peak around r ~ 0.5, which is the value one expects if particles have to 
diffuse through an angle A9 ~ l/Fg upstream, and an angle A^ ~ 1 downstream between 
shock encounters. 
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times as measured in up/downstream rest frame 
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Figure 7. The distribution of up- and downstream residence times for particles interacting with the shock in the case of 
isotropic scattering in both the up- and downstream rest frames. This shows the distribution for particles with E > 30 X -Einji 
which ensures that any influence of the initial conditions has decayed. 

Table 3. The spectral slope s, the mean ratio of initial and final energy < Ef/E\ > in a crossing cycle upstream > downstream 

— > upstream, and the return probability Prot for the case of pure deflection upstream and strong scattering downstream of 
the shock. 



< Ei/Ei > 



Pr< 



10 

103 



2.28 ± 0.01 1.67 ± 0.01 0.516 ± 0.028 

2.30 ± 0.01 1.65 ± 0.01 0.519 ± 0.005 

2.31 ± 0.01 1.61 ± 0.01 0.533 ± 0.005 



6.2 Case of regular deflection upstream 

We now turn to the case where upstream particles are deflected by a uniform regular fleld 
B±. Table 3 gives the results of our simulations. 

The energy gain per crossing is less than in the case of pure scattering, (Ef/Ei) ^ 1.64, but 
the return probability is higher, V^et ~ 0.523. As a result, the spectral index of the energy 
distribution of accelerated particles is only slighly steeper than in the case of pure scattering: 
s ^ 2.3. The results are internally consistent: from the average values of Table 2 we find 
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Figure 8. The angular distribution of particles in the downstream rest frame as a function of cos 6^ = /I. The smooth curve is 
the distribution of the pure scattering case derived from the results of Kirk et al. (2000), while the histogram is the simulated 
distribution in the case of regular deflection upstream and pure scattering downstream. The edge of the loss cone is located at 
/I = §. In the case of regular deflection upstream, the distribution extends only to JI > — i in agreement with the calculations 
presented in Section 4, Eqn. (bq). This curve is for particles with an energy E > 30 X -Binj, where the initial conditions no 
longer influence the distribution. 



1 + 



ln(l/P, 



ret; 



2.31 ±0.02 



In {Ei/E,) 

The differences with the case of pure scattering can be explained by looking at the angular 
distributions. Figure 8 gives the downstream angular distribution for pure scattering up- and 
downstream, and for the case of deflection upstream and scattering downstream. 

The angular distribution in the case of upstream deflection by a regular field is narrower, 
and cuts off at /i = — | (see Eqn. ^). This has two effects: 

(i) In the case of regular deflection, the downstream turning angle A6 = \6^d — ^-,u| 
between shock encounters is smaller, leading to a smaller value of (Ef/Ei) in an upstream- 
downstream- upstream crossing cycle; 

(ii) Since in this case there is (on average) a smaller fraction of particles with a large 
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Figure 9. The distribution of the energy gain, g = E{/Ei, in an upstream ^=> downstream =^ upstream crossing cycle. The 
curves shown are for particles with B > 30 X E^^j where the initial conditions no longer influence the dustribution. The thick 
line corresponds to the case of pure isotropic scattering on both sides of the shock, while the thin line corresponds to the case 
where particles are deflected by a regular upstream field, and are isotropically scattered downstream. 

value of ^^d, the chance Pret of re-crossing the shock into the upstream medium is larger, 
and the escape probability Pcsc correspondingly smaller. 

This explains the difference of the values in Table 2 and Table 3. If one plots the energy 
gain per crossing in a upstream-downstream-upstream cycle for both cases (figure 9) the 
difference is obvious: the distribution extends to higher values of Ei/ E\ in the case of pure 
scattering. 

Figure 10 shows the distribution of the upstream residence time tu{E) and the down- 
stream residence time t^ for Fg = 1000 and upstream defiection, as a function of the dimen- 
sionless variable 
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ZeB±TsCty^ 



upstream 

E 



T = < 



Vgid downstream 

The upstream residence time cuts off sharply at r = 3, which corresponds to 
3E 3 



t„ 



ZeB^cTs fi^Fs 



the behaviour found in our analytical calculations of Section 4 (see also Figure 2). The 

upstream residence time distribution peaks around 
E 

tu 



ZeB^cTs 
while downstream there is a rather broad distribution centered around r ^ 0.2, which 

corresponds to 

td~0.2 (VeY^ , 

a smaller value for t^ than obtained in the case of pure up- and downstream scattering. This 
latter result again shows that the average value of A^ = \6^d — 0^u\ is smaller in the case 
of regular deflection upstream. 

6.3 Initial shock encounter 

We have also simulated the initial shock encounter in order to calculate the energy gain at 
this first encounter, as well as the return probability Pjnj in this case. The latter quantity is 
essentially the fraction of particles picked up for further acceleration, and can therefore be 
considered as an injection probability for true (Fermi-type) relativistic shock acceleration. 
Note that our simulations assume that the shock itself is infinitely thin which, for these 
particles, may not be a very good approximation. We also assume that the particles are 
scattered isotropically in the downstream flow. We have considered two cases: 

• 'Cold' upstream particles with Ei ^ mc^ and \v\ ^ c; 

• Relativistic particles with 7^1 distributed isotropically in momentum. 

The results in both cases are virtually indistinguishable within the accuracy allowed by the 
simulations. We find that 
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Figure 10. The distribution of up- and downstream residence times for particles interacting with the shock through deflection 
by a regular magnetic field upstream, and isotropic scattering downstream. 
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Roughly 10% of the particles crossing the shock for the first time return upstream. A similar 
result has been found by Bednarz (2000) using a somewhat different simulation method. 

This result is of particular importance for models of Gamma Ray Burst afterglows. 
There it is often assumed that all particles crossing the shock are accelerated to a power- 
law distribution. Our results show that this assumption is probably too optimistic unless 
the thermalization of the bulk of the matter in the shock front also leads to a power-law 
distribution in energy for the downstream particles. 



7 CONCLUSIONS 

In this paper we have considered the process of shock acceleration at ultra-relativistic shocks 
with a shock speed Vs = /SgC such that the corresponding Lorentz factor satisfies Fg ^ 1. 

We have shown using both analytical methods and simulations that true Fermi-type 
shock acceleration, where particles cycle across the shock repeatedly, is subject to the effects 
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of relativistic beaming. Upstream particles are confined to an angle 9 < 2/rs with respect to 
the shock normal. This limits the energy gain per cycle so that the average ratio of final and 
initial energies in a cycle is Ei/E\ ~ 1.6 — 2.0, with about 50-60% of the particles escaping 
downstream each cycle. 

Like its non- relativistic counterpart, the mechanism is suprisingly robust in the sense 
that the energy distribution of the accelerated particles depends only weakly on the precise 
mechanism that confines particles near the shock. Downstream confinement always requires 
strong scattering in the momentum direction in order to allow particles to return to the 
shock. Upstream confinement can be achieved in two ways: by isotropic scattering of the 
momentum direction, or by deflection by a (quasi-) regular magnetic field. We find that in the 
first case the slope of the distribution equals s ~ 2.23, while in the latter case is is somewhat 
steeper, s ^ 2.30. Our results in the case of up- and downstream isotropic scattering agree 
with the analytical results of Kirk et al. (2000), and, in the limit Fg ^ 1, with the simulations 
of Bednarz & Ostrowski (1998) who employ a somewhat different simulation method. 

The largest particle energies can be achieved in the regime where particles are deflected 
rather than scattered upstream. This deflection can be due to a homogeneous upstream field 
or due to the sloshing motion in low-frequency MHD waves. The maximum energy attainable 
for a particle with charge Ze at a spherical shock of size R^ in the absence of losses is 

^^ax ~ ZeB TA , 

where B is the relevant upstream field strength, which is the component -B_l of the field 
perpendicular to the shock normal for a homogeneous upstream field, or the magnetic am- 
plitude 6B of the waves in the case that the upstream field is closely aligned with the shock 
normal so that Bj_ ^ B^/Ts, provided this amplitude is sufficiently large: 6B>B\\/rs. 

The momentum direction upstream will be able to diffuse between shock crossings if the 
upstream field is random on a scale ic < Rs for particles with an energy E^, > ZeB Tgic- In 
this case the maximum attainable energy, again in the absence of losses, is 



-Emax ~ ZeB Tg \J Rdc , 

with B the typical amplitude of the random field. 

The above two limits apply mostly to hadrons which suffer little or no radiation losses. 
For electrons (or positrons) synchrotron losses can limit the maximum energy to an energy 
of order 
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\o"t4b -d±/ 
where ^b is the downstream field amphfication factor defined in Eqn. (|iD|). All these energies 

are measured in the URF. 

The critical frequency z/g of the synchrotron radiation generated by electrons with Lorentz 
factor 7e = E/mec"^ is (Rybicki & Lightman, 1979, Ch. 6) 
eBi\ 2 



Att KnieC/ 

If synchrotron losses determine -Emax, particles with the maximum energy radiate photons 
with energy 

^^max ^ ^ (^] ^ 25 r, MeV , 
Stt \ a J 

in the URF, with a = e^ /he ~ 1/137.04 the fine-structure constant. This limit (roughly) 

applies to synchrotron photons originating in the upstream fiow, and to the Doppler-boosted 

photons originating in the downstream fiow. Note that the maximum photon energy is 

independent of the magnetic field strength in the source. 

We also considered the effect of the shock on particles at the first encounter. Assuming 
an infinitely thin shock with isotropic scattering downstream, we find that about 10% of the 
particles manage to return upstream so that they can participate in the shock acceleration 
process. They return upstream with an energy E ~ {].9T1Eq, where Eq is their energy before 
the encounter with the shock. This provides a natural injection process for shock acceleration 
at ultra-relativistic shocks. 

We note here that electrons will have to be pre-accelerated in order to be picked up by 
the shock acceleration process, assuming the thickness of the (collisionless) shock is of order 
the gyration radius of the downstream protons, which have an energy 

in the downstream rest frame. Cold upstream electrons entering the downstream fiow have 
an energy E^. ~ FgrneC^, and must be pre-accelerated by a different mechanism up to proton 
energies before they will start to see the shock as a discontinuity. A possible mechanism 
which could provide the necessary pre-acceleration has been proposed in the context of 
relativistic shocks in a pulsar wind by Hoshino et al. (1992). 
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